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FOREWORD 


This  work  was  performed  at  the  Supersonic  Wind-Tunnel  Laboratory  of 
Montana  State  University  by  A.  Demetriades,  Professor  Mechanical 
Engineering,  under  Grant  AFOSR  80-0267  with  the  Air  Force  Office  of 
Scientific  Research  of  Washington  D.C.  It  Is  part  of  a  series  of 
theoretical  and  experimental  studies  of  the  high-speed  mixing  character¬ 
istics  of  laminar,  transitional,  and  turbulent  shear  flows. 

In  a  practical  sense,  this  study  furnishes  theoretical  predictions 
of  wake  flows  needed  to  optimize  gas-dynamic  and  chemical  laser  design. 

As  such,  it  was  most  recently  motivated  by  work  done  by  the  author  for 
the  Air  Force  Weapons  Laboratory  and  TETRA  Corporation  of  Albuquerque, 

New  Mexico.  Specifically,  the  present  problem  was  first  discussed  in 
TETRA  report  TR-81-005  by  the  present  author,  in  which  the  solutions  for 
the  velocity  were  first  shown.  In  the  present  report,  the  subject  is 
continued  by  presenting  the  solutions  for  the  thermodynamic  variables  as 
well.  For  the  sake  of  completeness,  the  earlier  velocity-field  solutions 
are  also  briefly  discussed. 

A  second,  related  motivation  arises  from  the  need  to  prepare  theoretical 
predictions  for  certain  experiments  planned  in  the  laboratory's  supersonic 
tunnel  for  the  AFOSR  program.  These  experiments  feature  the  basic  fluid 
mechanics  of  asymmetric  wakes  and  shear  layers  at  supersonic  speeds,  for 
which  predictive  formulas  are  presented  here.  Clearly,  if  the  latter  are 
verified  by  the  forthcoming  experiments,  the  present  analysis  can  become 
an  important  tool  in  laser  cavity  design. 

Encouragement  of,  and  interest  in  this  work,  has  therefore  come 
from  diverse  sources  which  the  author  wishes  to  acknowledge  here.  The 
assistance  and  cooperation  of  Drs.  W,  Moeny  of  TETRA  and  P.  J.  Ortwerth  and 
L.  Wilson  of  AFWL,  as  well  as  of  Capt.  M.  Francis  of  AFOSR  is  noted  with 
gratitude. 
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1 .  Introduction 


The  problem  addressed  In  this  report  Is  the  steady,  laminar  two- 
dimensional  wake  formed  by  a  homogeneous  fluid  streaming  past  an  Infinitely 
thin  plate.  As  enunciated,  this  problem  was  first  solved  for  Incompressible 
flow  ty  Tollmlen  (Reference  1)  and  Goldstein  (Reference  2).  These  solutions 
take  a  simple  self-similar  form  far  away  from  the  trailing  edge  (T.E.),  which 
has  also  been  used  as  a  convenient  starting  point  in  the  discussion  of 
compressible  wakes  (Reference  3).  In  the  region  Immediately  past  the 
T.E.  the  classic  solutions  of  References  1  and  2,  however,  have  been 
arrived  at  in  the  form  of  matching  expansions,  which  produces  tabulated 
numerical  results  rather  than  a  single  analytic  expression.  It  would 
clearly  be  desireable  to  have  closed-form  analytic  solutions  instead,  In 
a  way  Improving  the  perception  of  the  wake  behavior. 

This  report  sets  out  to  accomplish  two  objectives,  the  first  of 
which  is  the  derivation  of  formulas  for  the  wake  properties  which  are 
valid  from  the  T.E.  to  an  infinite  distance  downstream.  This  Is  done 
at  any  stream  Mach  and  Reynolds  numbers  and  any  T.E.  surface  temperature 
relative  to  the  stagnation  temperature,  by  restricting  the  analysis  to 
isobaric  flow,  Prandtl  number  one  and  constant  Chapman-Rubesln  factor. 

These  restrictions  are  certainly  not  severe  enough  to  obscure  the  dynamics, 
and  can  be  removed  once  the  basic  problem  is  understood. 

The  second  objective  is  to  extend  the  solutions  to  a  wake  with 
initial  asymmetry.  In  reallistic  terms,  such  an  asymmetry  is  more  likely 
to  occur,  say,  when  two  adjacent  nozzles  discharge  parallel  streams,  as 
shown  on  Figure  1,  than  in  the  classic  flat-plate  wake.  In  fact,  the 
motivation  for  this  work  arose  from  questions  of  gas-dynamic  and  chemical 
laser  design  optimization,  which  deal  with  geometries  such  as  that  of 
Figure  1.  Since  the  two  streams  are  Identical,  the  asymmetry  can  arise 
only  due  to  surface  conditions  at  the  T.E.,  such  as  different  momentum 
thicknesses  or  different  wall  temperatures  on  either  side  of  the  T.E. 

The  latter  is  rather  unreal  11  Stic,  however,  so  long  as  we  consider  the 
partition  near  the  T.E  as  "infinitely  thin".  Thus,  we  consider  only  the 
asymmetry  due  to  momentum- thickness  differences  in  the  two  merging 
boundary  layers. 
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It  must  be  '.toted  that  the  method  of  solution  for  both  the  symmetric 
and  asymmetric  cases  rests  on  a  linearization  assumption.  However,  since 
the  initial  and  asymptotic  conditions  at  large  distances  are  recovered, 
the  penalty  for  this  assumption  Is  nowhere  evident. 

2.  Statement  of  the  Problem  and  Definition  of  Terms 


The  solution  of  the  wake  problem  given  herein  is  subject  to  the 
following  restrictions: 

1)  The  flow  is  everywhere  parallel  and  two-dimensional. 

2)  The  pressure  is  everywhere  constant. 

3)  The  fluid  is  chemically  homogeneous  and  non-reacting. 

4)  No  volume  effects  (e.g.  volume  heat  addition  or  body  force)  are 
allowed. 

5)  The  flow  is  laminar  and  steady. 

A  lesser  assumption  is  that  the  Prandtl  number  is  unity  everywhere, 
chosen  consistently  with  the  use  of  the  Crocco  relation  for  the  temperature 
distribution  in  the  boundary  layer  at  the  trailing  edge  (T.E.). 

Within  these  restrictions,  the  solutions  are  provided  for: 

(a)  Any  stream  Mach  number  M  . 

(b)  Any  Reynolds  number  Re. 

(c)  Any  wall-to-stagnation  temperature  ratio  T  /T  (T  is  the  T.E. 

w  u  w 

surface  temperature) , 

(d)  Any  asymmetry  ratio  ©i /©2  =  P. 

The  solutions  are  thus  provided  as  functions  of  x,  y,  M,  T  /T  ,  P 

W  0 

and  Re,  the  latter  implicitly  through  the  definition  of  the  non-dimensional 
variables.  Wh.  =  Gg  or  P  =  1 ,  the  solution  is  termed  symmetric,  whereas 
the  term  "completely  asymmetric"  denotes  the  solution  for  P  =  0  (or  P  *  »). 
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The  terms  defect,  deficit,  and  decrement  are  used  to  denote  characteris- 


A 


s 

I 


► 


f 

* 


-  r 


c 


tic  groupings  of  the  flow  variables  at  a  certain  x  from  the  T.E.  For  example, 
for  the  temperature: 


T(0) 

-  Te  or  'c  -  Te  : 

deficit 

T(0)  -  Te 

or 

Tc  *  Te  : 

defect 

Te 

T(0)  -  Te 

or 

Te 

Tc  -  Te  : 

decrement 

VTe 

VTe 

The  subscript  "c"  (for  "center")  is  used  to  denote  the  extremum  of  the 
temperature  profile  at  a  certain  x  for  the  general  case  of  wake  asymmetry 
(P  ^  1).  When  P  =  1 ,  Tc  coincides  with  T(0),  the  latter  symbolism  meaning 
the  temperature  at  y  =  0  which  is  also  the  wake  "center"  (plane  of  symmetry) 
when  P  =  1 . 

The  term  "profile"  is  reserved  for  specific  non-dimensional  groupings 
such  as 


ue  -  u(0)  ue  -  uc 

while  the  term  "variation"  is  often  used  to  represent  groupings  such  as 

1  -  U/Ue. 

3.  Solution  of  the  Problem  for  Arbritrary  Asymmetry 
3.1  The  Method  of  Solution 


The  present  objective  is  to  find  the  fluid  properties  (velocity, 
temperature,  density,  etc.)  in  the  wake  as  a  function  of  the  coordinates, 
the  flow  parameters,  and  the  asymmetry  ratio  P.  To  this  end,  we  utilize 
the  solution  reported  by  Gold  (Reference  4  )  for  a  W3ke  beginning  at  the 
T.E.  with  profiles 


‘(<,1  )  3 

3 


u 

u* 


(1) 


(2) 


Ar 


h)  = 


T  (A**) 


Te 


of  the  velocity  and  temperature.  Gold  employed  the  Oseen  approximation 
to  the  boundary -layer  equations  to  find  the  solutions  as  follows: 

=  ■£7='  |  uo  ^ l  o) 

_  cc— 

jL=  j  (4) 

-00 


where  uQ,  h0-'  initial  velocity  and  temperature  profiles  ,  and 


X  s 

X* 

(5) 

L 

**  - 

±{r£l 

L 

(6) 

If  * 

u  * 

!  L  <*•» * 

Jo  h  1 

(7) 

L  =  characteristic  length,  to  be  derived. 

3  U-%L  (8) 

G-fi > 31  7^  'exP  (9) 


The  Oseen  approximation  is  a  linearization  restricted  to: 


U  <  <  I 
•fc.  <  <  i 


(10) 

(ID 


Note  that  dimensional  properties  above,  are  usually  denoted  by  an 
asterisk  (*)  or  a  subscript  "e". 
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Three  additional  remarks  must  be  made  as  regards  the  above  solution. 
First,  Gold  assumed  the  Chapman-Rubesin  factor  as  constant;  second,  the 
coordinate  x  appearing  in  eq.  (5)  was  written  by  Gold  as  x/r,  to  account 
for  non-unity  Prandtl  number.  Here  we  set  S'  =  1  to  be  consistent  with  our 
later  use  of  the  Crocco  relation.  This  difference  can  be  easily  accounted 
for  in  numerical  calculations,  and  is  in  anyway  a  rather  unimportant  feature 
of  this  analysis. 

Third,  Gold  calculated  that  the  solutions  (3)  and  (4)  are  rather 
insensitive  to  the  initial  profiles  and  FQ  especially  far  from  the 
origin,  and  that  the  governing  factor  is  the  initial  momentum  thickness 
rather  than  the  initial  profile  shape.  This  should  be  kept  in  mind  in 
the  next  section. 

3.2  The  Initial  Profiles 


If  the  shed  boundary  layers  have  momentum  thicknesses  6.i  and  02 
respectively,  then  the  initial  velocity  profile  at  the  T.E.  for  purposes  of 
calculating  the  wake  development  is: 


For  j  >  0 


a.  -*•* 


ue 


^  >  1 


with  k-j ,  k^  evaluated  from  the  linearized  version 

o°_ 

e  ^  j  i  M  * 


°  3 


of  the  momentum  integral 


where 


~  3  P  ,  v- 

1  =  ?.•*» 


Use  of  (14)  Instead  of  (15)  as  a  boundary  condition  Is  consistent  with 
Gold's  theory  which  is  used  past  the  T.E.  and  it  implies  that 


* 

HI 

1 

C 

o 

A 

A 

(17) 

ue 

Using  (12), 

(13)  in  (14)  gives 

^1  =  91  ’  **2  62 

(18) 

so  that  the 

initial  velocity  profile  reads 

-Va, 

UQ  s  fi  =.« 

(19) 

-  -  .  J*'*'  _  cp* 

<  0  ^0  *  ^  ^ 

(20) 

with 

S 

(21) 

•p  *  e, 

(22) 

The  choice  for  an  exponential  initial  velocity  oroflle  is  arbritrary, 
although  it  has  obvious  qualitative  similarities  to  a  "typical"  profile  as 
Figure  2  shows;  clearly,  too,  it  is  algebraically  simple  and  should 
therefore  expedite  the  solution.  The  earlier  comments  should  also  be 
recalled,  regarding  the  insignificance  of  the  initial  profile  details  to 
the  wake  development.  It  is  more  important  to  note  that  the  chosen 
profiles  (19)  and  (20)  cannot  possibly  fulfill  the  condition  (10).  We 
shall  see  from  the  results  that  the  solution  is  analytic  for  all  x  >  0, 
however,  and  that  a  very  reasonable  continuous  solution  is  obtained. 

The  parameter  P,  called  the  asymmetry  ratio,  is  a  key  parameter  in 
this  analysis  because  it  shows  the  initial  "skewness"  of  the  wake. 

Because  of  the  interchangeability  of  01  and  e2  the  range  0  <  P  <  1  covers 
all  possible  ratios  of  and  e2< 
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I 


For  consistency,  the  definition  of  y  must  be  carried  from  the  T.E. 
to  all  downstream  positions.  This  allows  us  to  define  the  length  scale 
L  by  combining  (6)  and  (21): 


1  . 

1 

e, 

(23) 

« 

(24) 

I 


so  that  the  non-dimensional  variables  x  and  y  in  the  proposed  solution 
(3)  become,  according  to  (23)  and  (24): 

v  =  -  (25) 


* 

9. 


(26) 


We  perceive  that  the  choice  of  L  can  be  further  improved,  to  avoid 
awkward  mathematics  such  as  for  F  =  0  (6^  =0).  Besides,  it  should  be 
clear  that  the  wake  should  be  controlled  by  Its  total  drag,  i.e.  0-j  and 
92,  and  by  P  which  shows  how  this  drag  was  initially  distributed  between 
streams  1  and  2.  We  will  make  this  Improvement  on  L  later,  since  it  Is 
more  convenient  to  first  present  the  solution  in  terms  of  the  coordinates 
as  defined  by  (25)  and  (26). 


The  initial  temperature  profiles  were 

found  with 

the  aid  of  the 

Crocco 

relation: 

1  i 

~r  tv  /  y-t  Tw  \  U 

±  +(  H-V-  M*  ~  ^  )  -  - 

Tt  Te  2-  T«  Me 

V-i . .  i 
a 

U  *■ 

tW« 

(27) 

f 

Using  the  definitions 

(28) 

B  S 

(29) 

» 

C  e  -  ~ 

(30) 

7 


(31) 


as  well  as  equations  (19)  and  (20)  we  obtain 

>}>«>■•  £0  =  C 

-  T-Pu 

Co:  k  =  B  -e  ^  4  C  € 

o 

The  Crocco  relation  {21)  is  valid  when  the  Prandtl  number  €-  1.  For 
this  reason,  as  already  mentioned,  consistency  requires  that  the  wake 
solutions  are  found  for  €  =  1  also,  as  we  intend  to  do. 

Typical  initial  temperature  profiles  according  to  (31)  and  (32),  are 

illustrated  on  Figure  3,  and  are  useful  in  stressing  two  points.  First, 

the  examples  of  Figure  3  were  chosen  to  illustrate  the  possibility  that 

for  certain  M  and  T  ,/T  we  can  obtain  A  >  >  1,  which  violates  the  assump- 
e  w  o  o 

tion  (11).  This  violation  is  more  serious  than  in  the  c?se  of  the  velocity  I 

where  uQ  s  1  at  most;  thus  the  temperature  results  of  the  present  analysis 

should  be  viewed  with  caution  near  x  =  0.  Second,  the  distribution  of 

temperature  is  plotted  in  two  different  ways,  the  "profile"  function  being 

shown  on  vho  right  of  the  Figure.  The  algebraic  form  of  the  profile  often  $ 

creates  unusual  curves,  a  point  to  be  kept  in  mind  during  the  presentation 

of  results. 


i' 

(32) 


3.3  Solutions 

If  the  initial  profiles  (19),  (20),  (31)  and  (32)  are  substituted 
in  eq.  (30)  and  (4),  simple,  closed-form  solutions  can  be  Immediately 
obtained  for  i7  and  : 

u(y^  jP)  =  4  >)  ¥ 

*  <33) 

s[’B[e,‘'*(i-64reo-2i»+e*  *  Vi-  (34) 

-E-t^Vx  (^♦'s))]  +  cIe<,-z»(i-E^rt^-i 

As  expected,  u"  does  not  explicitly  depend  on  Mg  and  TW/TQ,  although  it 
does  so  implicitly  through  the  transformed  variable 
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The  solutions  (33)  and  (34)  have  been  successfully  subjected  to  the 
following  tests: 

(a)  At  x  *  0,  the  Initial  profiles  (19)  etc.  are  recovered. 

(b)  At  y  =  ± 03  ,  u  *  as  expected. 

(c)  At  x  a®,  the  familiar  Tollmlen-Goldsteln  asymptotic  forms  are 
recovered  as  expected,  Independent  of  P.  (References  1  and  2) 
The  Tollmlen-Goldsteln  asymptotes  will  be  presented  later. 

(d)  If  01  and  02  are  interchanged  relative  to  the  positive  sense  of 
y,  then  the  asymmetric  solutions  (PjM)  also  interchange,  as  they 
should. 


Solutions  (33)  and  (34)  are  sufficient  to  determine  the  entire  flow 
field  since  the  pressure  is  constant;  thus  the  density,  Mach  number,  etc. 
can  be  determined  directly. 

Before  applying  the  solutions  to  special  cases.  It  is  here  appropriate 
to  change  the  independent  variables  into  a  more  meaningful  form.  The 
total  momentum  thickness  of  the  wake  is 

(S)  =  9(  +  0*.  (35) 


and  it  appears  suitable  to  use  ®  in  place  of  61  in  the  definitions  of  x 
and  y  (see  (25)  and  (26)).  The  new  variables  therefore,  are: 


& 


V*  - 

o  © 

The  connection  between  x,  y,  and  x'. 


CP+O' 


x* 

es)rf?e'  (36) 

(37) 

'  can  be  found  using  (35): 

(38) 

(39) 
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The  solutions  (33)  and  (34)  thus  become: 

A.# 


and 


+  eO')V-^,-  (  +  £^(^.Pf/.o)] 


(40) 


;  (41) 


•*>  *  4{SI>  - - T'  ^  *«V-  Vlli'f-zf&V 

((P(lW,+  y^+c[e  r 

d-s-,^'?1  -2^,V)  +e*^"'x'',’l<,‘")^-s4(2(t>*.)r,'+^]  j 

It  should  be  parenthetically  noted  that  the  variable  is  insensitive 

to  the  choice  of  the  normalizing  length,  as  can  be  seen  by  inspection  of 
eqs.  (36)  and  (37). 

4.  The  Symmetric  Wake 

As  mentioned  in  the  Introduction,  a  complete  analytic  solution  from  the 
T.E.  downstream  has  not  yet  been  reported  for  the  symmetric  two-dimensional 
wake,  regardless  of  M  or  TW/TQ.  Such  a  solution  can  be  obtained  directly 
from  (40)  and  (41)  by  setting  P  *  1.  Then, 

-  - - «4» -**’<.♦«* it,  _,<!.»]  w 

;  =  i[B(44y,-^(l-^c^x--^)  +  e4AV<(.E,+  u;w^,i)) 

f  c  (<(*'-2^))  +  e'4x  ++'l(|'64^'t^,)1)](43) 


These  solutions  are  naturally  symmetric  in  y'. 


2/*' 


In  the  symmetric  case,  the  wake  center  occurs  at  y'  3  0,  (i.e.  the  plane 
of  symmetry)  so  that  if  the  wake  defect  is  defined  as 

ue-U<o)  u-uc 


W  ~ 


U, 


=  w  (*') 


(44) 


and  the  velocity  profile  as  usual  by 

rJ 


r*t  Uo  -U  ~  >  ,  \ 

U  =  - -  -  U  ^  ) 

Ue-u^o) 


(45) 
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\ 


then 


Ui  =■  WU 


(46) 


I 


( 


t 


t 


It  can  be  easily  shown  from  (42)  that  the  defect  Is 

=  i)  -  e  0-En£2 (47) 

and  that  W(x"  s  0;  P  s  1)  a  1  as  expected.  At  very  large  this  defect 
becomes 

w(y'-“  if*°  =  jfe*  <48> 

which  is  well-known  of  ths  Tollmlen  (Ref.  1  )  and  Goldstein  (Ref.  ,*>  )  and  which 
has  since  been  discussed  sy  numerous  authors  (Ref.  3,  4). 


The  velocity  proiile  begins,  at  x"  ■  0,  with  the  double-cusp  shape  given 
by  (19)  and  (20)  and  P  «  1.  The  discontinuity  at  y'  a  0  disappears  just  as 
soon  as  the  fluid  leaves  the  T.E.;  for  all  x  >  0  the  u,  found  by  combining 
(42),  (46),  and  (47),  is  continuous  and  analytic  In  y'.  At  largo  x'  it 

•x. 


can  be  stiown  from  (45)  that 

U  Cv'-oo  3?»i)  -  e 


(49) 


This  simple  result  has  been  long  known  as  the  asymptotic  form  of  the 

profile,  which  depends  on  the  similarity  variable 

/ 

Zfy.' 


F  = 


(50) 


'he  temperature  distribution  at  the  T.E.  as  given  by  (31)  and  (32)  is 
of  course  generally  discontinuous  at  y  *  y'  s  0.  Just  like  the  velocity, 
this  distribution  develops  a  "rounded  top"  and  Is  everywhere  analytic 
Immediately  past  the  T.E.  according  to  (43).  On  the  plane  of  symmetry 
y'  *  0,  the  temperature  defect 

T  (q) 


±  s  <(B)  »  sP*  I  j  !*€,■&)  -  * 

*0  ® 

=  Be4*0  +-  c  4/xO 

At  the  trailing  edge  (x-  *  0)  the  Initial  defect  is 

o)  =  3+C 

Te 

11 


(51) 


(52) 


as  also  Indicated  by  (31)  and  (32).  Far  from  the  T.E.  on  the  other  hand, 

i  (V ^  (53) 

Combining  with  (47)  we  note  that  at  large  x'  the  ratio  of  the  temperature 
to  the  velocity  defects  is  a  constant  dependent  on  and  TW/TQ  only: 

±=B+£-rf.M:(.+^)+r„_i}  (M) 

Here  the  first  term  on  the  r.h.s.  Indicates  the  compressibility  effect  on 
wake  temperature,  while  the  second  term  Indicates  the  effect  of  the  surface 
temperature  only.  There  Is,  of  course,  a  Prandtl  number  effect  which  is 
Ignored  hpre  as  per  earlier  comments;  If  the  Frandtl  number  «*  was  Included, 
then  the  entire  r.h.s.  of  (54)  would  be  multiplied  by  C  . 

In  addition  to  the  temperature  defect  as  defined  In  (51),  we  can 
consider  the  temperature  decrement 

»_  T(o)-Te  ..  Tv.v 

*  TWi  '  *  j  (55) 

which  is  always  unity  at  the  T.E.  as  opposed  to  the  defect  which  can  be 
larger  than  unity  at  x"  =  0.  The  decrement,  which  is  simply  related  to  t 
via  the  ratio  Tw/TQ,  is  useful  when  the  latter  is  different  from  unity. 


The  temperature  profile 


T  5 


T-Te 

T(o)-Te 


is  equal  to  the  ratio  ^  /t  and  can  be  thus  found  from  (43)  and  (51).  At 
the  T.E.,  the  profile  is  given  for  P  *  1  by  the  Initial  conditions  (31)  and 
(32);  at  large  x%  on  the  other  hand,  the  profile  asymptotic  form  Is  again 


T(Aco)  -  € 


=  e  “  U  -OD) 


which  is  again  a  consequence  of  1.  Parenthetically  note  that,  In  general. 


(y-cD  -  £ 


-CF' 
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(58) 


I 


I 


< 


The  fact  that  the  asymptotic  velocity  and  temperature  profiles  (49)  and 
(57)  agree  with  the  long-known  results  of  Tollmlen  and  others  Is  an  Important 
cross-check  of  the  present  analysis. 

Other  flow  variables,  such  as  the  density,  total  temperature  and 
Mach  number,  can  be  obtained  In  the  above  manner  quite  directly.  For 
the  sake  of  completeness,  these  variables  will  be  given  below. 


Using  Its  customary  definition,  the  density  defect  Is 

*  s  UJL*  m  _±_ 

fe  +  ' 

with  t  given  by  (51).  Similarly,  the  density  decrement  Is 


The  density  profilers 

h-f 


a 


fc,  +  1 

T  »- - 

Ti  +•  1 


(59) 


(60) 

(61) 


An  Inspection  of  these  formulas  for  the  density  reveals  once  again  that 
the  density  "lags"  the  fluid  temperature  In  the  attainment  of  similarity  at 
large  x'. 


The  local  Mach  number,  M,  In  the  wake  is  given  by 

l-uw 


M  35  MC 


TiTft 


and  on  the  plane  of  synnetry  ( u  -7  -  I  ): 

/  x  1 

M  (0)  »  t =r 

e  VTPt 


(62) 


(63) 


Finally,  the  total  temperature  defect  ^ 

/  +0.2  M(o)  . 

-I  (r^l‘4) 


T0(o) 

Tee 


H  =  =  l  *■+») 


1  +  0.2  Me 


while  the  total  temperature  decrement  is 

-  _  TM-T.e  _  H 

n  =  ’  1 


TW 


-  I 
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(64) 


Tw-7oe 


(65) 


c  r-i. 


and  the  corresponding  profile  Is 

~  B  T.  -Tot  =  a*?*) 

o)-7ot  (  i+±) 


I40.2M 
I  +  0.2  Mfr  ' 
i  ^o.iwa)1 
\  i-o.imV 


(66) 


Note  that  In  the  last  few  formulas,  above,  dealing  with  the  local  total 
temperature,  the  symbol  for  the  latter  Is  TQ.  The  quantity  TQe  Is  the  glven 
total  temperature  of  the  stream,  which  In  earlier  sections  was  denoted  by 


Sample  computations  and  plots  for  the  symmetric  wake  will  be  presented 
in  Section  5. 


5.  The  Completely  Asymmetric  Wake 

Early  computations  made  In  this  work  with  eqs.  (40)  and  (41) 
showed  that  considerable  departures  from  the  symmetric  wake  case  given 
above,  occur  only  for  extreme  values  of  the  asymmetry  ratio  P,  that  Is  for 
9l  >  >  e2  or  e1  <  <  e^.  It  therefore  seems  reasonable  to  limit  the  analysis 
of  wake  asymmetry  to  the  extreme  case  6-j  *  0  (P  »  0)  or  02  *  0  (P  =>  °°). 

These  two  extreme  cases  are  of  course  Identical  since  one  can  be  obtained 
from  the  other  by  mirror  reflection  of  Its  flow  field  about  the  y'  *  0 
plane.  We  shall  therefore  discuss  only  the  P  *  0  case,  which  means  that 
one  of  the  two  merging  boundary  layers  is  vanishingly  small  relative  to 
the  other.  We  call  this  the  "completely  asymmetric  wake".  (Of  course, 
the  general  solutions  (40)  etc.  can  be  used  for  arbritrary  asymmetry). 


For  P  *  0,  the 

initial  profiles  at  the 

T.E.  are: 

^>°  : 

I 

~  0  ~ 

5/<0  / 

(67) 

uj.  <  0  : 

lA0  ^  € 

€  *  e 

(68) 

V>°  : 

k0  -  0 

(69) 
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(70) 


<  0 


i.„  =  Be  a'+  c  e 


cel>' 


It  Is  both  necessary  and  Intuitively  possible,  here,  to  preview  the 
development  of  this  asymmetric  wake  beyond  the  trailing  edge.  As  x' 
Increases,  the  above  Initial  profiles  will  tend  to  "round  off",  but  the 
"center  of  gravity"  of  the  profile  Initially  lies  In  the  y*  <  0  half  space. 
It  Is  therefore  likely  that  the  profile  maximum  (actual  velocity  minimum) 
will  lie  In  y'  <  0  and  not  y“  ■  0  as  was  the  case  of  wake  symmetry.  This 
maximum  will  be  referred  to  as  the  "wake  center",  and  Its  path  (the 
"deflection")  away  from  y'  •  0  Is  of  some  Interest.  For  '■  ample,  we  can 
Intuitively  predict  that  the  wake  symmetrizes  at  large  x',  but  not  whether 
the  center  will  return  to  the  y*  «  0  plane.  The  subscript  "c"  for  "center" 
will  In  anyway  assume  the  role  the  symbol  (0)  had  for  the  symmetric  case 
(see  Section  2  also). 


The  solutions  for  the  completely  asymmetric  wake  can  be  directly 
obtained  from  (40)  and  (41)  In  the  limit  P  «  0: 

V)P=o e  'VsHfl/x+I/)  +|e  (72) 


(71) 


These  equations  exhibit  the  proper  limiting  behavior;  for  example, 
in  the  limit  x'  ■  0,  eq.  (71)  reduces  to  (67)  and  (68),  while  in  the  limit 
x"  «  «,  It  reduces  to  (48).  It  is  equally  easy  to  show  that  equations  (71) 
and  (72)  are  asymmetric  relative  to  y",  and  that  at  y'<  a  0  the  iu  obtained  from 
(71 ) ,  that  is 


U  (x'0'xo')P=o)  =  4  «•* 


(73) 


does  not  represent  the  apex  of  the  profile  as  it  did  In  the  symmetric 
case  (In  References  5  and  6,  eq.  (73)  was  called  the  "pseudodefect"). 
In  the  present  case,  the  velocity  defect  is  defined  as 


W  =  I  -  ^  =  w  (/  o  ;  P-  0) 


(74) 
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where  uc  is  the  minimum  velocity  at  any  x',  herein  also  called  the  "center" 
velocity.  Then  the  velocity  profile  is 


(A 


U  e-Uc 


(75) 


The  center  velocity  uc  as  well  as  the  position  y*  (the  "center  position") 
at  which  it  occurs  cannot,  unfortunately,  be  derived  in  closed  form.  It 
will  therefore  be  discussed  in  connection  with  the  computations  of  Section  7. 


The  discussion  of  the  asymmetric  temperature  field  follows  quite  closely 

the  above  discussion  for  the  velocity.  The  definition  of  the  temperature 

defect,  decrement,  profile,  etc.  follows  the  rule  that  the  symbol  (0)  used 

for  the  symmetric  wake  (eqs.  (45),  (56),  etc.)  is  now  changed  to  a  subscript 

c.  Note,  however,  that  the  apex  of  the  temperature  profile  is  again 

impossible  tc  derive  analytically  in  closed  form.  We  therefore  have  no 

'a  priori"  evidence  that  the  wake  center  found  from  the  temperature  profile 

coincides  with  the  wake  center  y\  ,  discussed  earlier  for  the  velocity 

cu 

profile.  This  point  will  be  clarified  in  the  computations  presented 
in  Section  7. 


The  remaining  properties  of  the  completely  asymmetric  wake  such  as 
the  density  can  also  be  presented  in  the  form  of  defects,  profiles,  etc. 
defined  as  for  the  symmetric  wake  but  with  subscript  "c"  in  place  of  the 
symbol  (0).  Results  will  be  shown  below. 

6.  Computations  and  Graphic  Results  for  the  Symmetric  Wake 

The  formulas  derived  in  Section  *  for  the  symmetric  wake  (P  *  1)  have 
been  utilized  to  present  graphically  the  behavior  of  the  fluid  properties 
in  such  a  wake  from  the  T.E.  (x'  =  0)  to  far  distances  downstream.  These, 
which  are  shown  on  Figures  4  through  16  and  which  will  be  discussed 
briefly  in  this  Section,  were  obtained  by  numerical  computations  with  the 
programs  shown  in  the  Appendices. 

The  chief  interest  here  lies  in  the  way  the  flow  adjusts  from  a 
discontinuous  profile  at  the  T.E.  to  a  continuous,  analytic  profile 
immediately  past  the  T.E.  We  are  also  interested  in  the  speed  with 
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which  the  asymptotic  solutions  (48),  (49),  etc.  are  attained.  We  are  no 
less  interested  in  the  effect  of  stream  Mach  number  M  and  of  the  temperature 

C 

ratio  Tw/T0  on  these  adjustments.  In  this  respect,  it  should  be  recalled 
that  the  fluid  velocity  is,  in  the  compressible  transformation  of  eq.  (16), 
unaffected  by  Mfi  or  TW/TQ,  but  that  the  thermodynamic  variables  are 
strongly  affected  by  these  parameters. 


It  is  well  to  keep  a  physical  perspective  of  the  distances  involved 
downstream  of  the  T.E.,  and  in  this  sense  the  kinematic  boundary  layer 
thickness  £  is  a  more  meaningful  length  measure.  For  example,  suppose 
we  are  interested  in  the  flow  in  a  gas-dynamic  laser  cavity  10  cm. 
downstream  of  the  nozzle  array.  Suppose,  further,  that  the  nozzle  exit 
Mach  number  is  high-supersonic  with  a  laminar  &  =  0.5  mm.  In  this  case, 
the  ratioS/0  is  of  order  20,  while  Re0  can  be  taken  to  be  400  (i.e. 
near  the  upper  limit  for  laminar  flow).  Then,  from  eq.  (36): 


X  = 


X1 


/  0 


4  400 

20 


=  0.2.5- 


(76) 


In  this  example  we  see,  therefore,  that  the  range  x'  <  1  is  of  major 
practical  interest. 


Figures  4  and  5  show  the  velocity  defect  according  to  eq.  (47).  The 
experimental  points  on  Figure  5  will  be  discussed  later,  but  the  comparison 
with  other  theories  is  in  order  at  this  point.  Specifically,  the  good 
agreement  with  Goldstein's  theory  (Ref.  2  )  should  be  noted,  as  well  as 
the  difference  from  the  "asymptotic  solution",  i.e.  eq.  (48).  The  latter 
seems  to  be  adequate  In  giving  W(x')  beyond  x'  -  1,  but  expectedly  fails 
to  do  so  at  x'  <  1.  Figure  4  demonstrates,  by  the  proximity  of  the  solid 
to  the  dashed  curves  that  for  very  short  distances  past  the  T.E.  the  profiles 
are  still  almost  identical  with  the  initial  profiles  uQ,  TQ  marked  "theory", 
(the  solid  curve  for  u  has  a  rounded  apex  not  visible  on  this  graph).  The 
evolution  of  the  velocity  profile  u  is  shown  on  Figure  6.  The  use  of 
two  different  abscissas  in  this  Figure  is  necessary  to  show  now  the  profile 
approaches  the  similarity  solution  (49),  occurring  past  x'  =  1  as  noted. 
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There  are  infinite  combinations  possible  of  Mg  and  TW/TQ  which  can 
be  used  to  examine  the  behavior  of  the  defects,  decrements,  and  profiles 
of  the  thermodynamic  variables  according  to  eqs.  (51),  (55),  (56),  etc. 
of  Section  4.  Generally,  computer  runs  with  the  programs  shown  in  the 
Apprendices  were  restricted  to  Mg  =  0,  2,  4,  and  6,  and  TW/TQ  values  of 
0.2,  0.6,  1,  1.4,  and  1.8  (the  "cooled",  "adiabatic"  and  "heated"  wake). 

The  temperature-decrement  plots  of  Figures  7  and  8  indicate  the 
significant  finding  that  the  wake  behavior  subdivides  naturally  into 
two  categories.  The  first  is  the  high-Me,  highly  cooled  type  of  wake 
(Figure  7  )  which  is  also  the  one  of  practical  importance  in  laser 
applications.  Here  the  decrement  is  "normal  for  Mg  <  4,  but  between 
Mg  =  4  and  5,  the  decrement  switches  sign — in  some  measure  due  to  the 
definition  of  the  decrement.  For  all  other  cases,  however,  the  behavior 
is  typically  shown  on  Figure  8  ;  the  decrement  decreases  monotonically 
from  the  T.E.,  and  is  always  included  within  the  band  defined  by  each 
pair  (solid  and  dashed)  of  curves.  For  example,  for  adiabatic  or  heated 
wakes  (T  /T  >1)  and  at  x"  =  1,  the  decrement  always  lies  between  0.25 

W  0 

and  0.35  regardless  of  Mg. 

The  effect  of  compressibility  and  heat  transfer  on  the  temperature 
profile  is  shown  by  the  three  plots  of  Figures  9,  10,  and  11.  The 
effect  of  high  Mg  alone  or  T.E.  cooling  alone  does  not  create  any 
unusual  behavior  in  t  as  Figures  9  and  10  indicate;  the  adiabatic  hyper¬ 
sonic  wake  is  simply  thicker  than  its  low-speed  cooled  counterpart. 

When  high  Mg  is  combined  with  cooling,  however,  Figure  11  shows  that  the 
profile  has  large  off-axis  peaks  in  the  temperature,  inherited  from  the 
non-monotonic  initial  temperature  distribution  at  the  T.E.  (see  dashed 
curves  in  Figure  3  and  also  in  the  upper  portion  of  Figure  11). 

Tne  density  and  total  temperature  defects  are  shown  on  Figures  12 
and  13.  Note  that  for  Mg  =  5  and  TW/TQ  =0.2,  the  density  defect  increases 
quickly  from  0.16  to  0.52  a  short  distance  (x'  =  0.05)  downstream  of  the 
T.E.  and  thereafter  decreases  monotonically. 
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Figures  14  and  15  illustrate  the  differences  among  velocity, 
temperature,  and  density  in  the  speed  with  which  each  tends  to  achieve 
the  asymptotic  state.  The  high-supersonic,  cooled-wake  environment  has 
been  selected  for  these  examples.  Figure  14  plots  the  three  defects 
together  with  their  asymptotes  (i.e.  eqs.  (48)  and  (53)),  since  the 
latter  are  sometimes  used  as  approximations  to  the  actual  defects.  We 
see  that  by  x'  =  0  the  actual  velocity  defects  are  within  15  percent 
of  their  asymptotic  approximations,  but  that  the  density  defect  is  closer 
to  one-half  of  what  its  asymptotic  formula  predicts.  Similarly,  for 
x'  -  0.1  Figure  15  shows  that  the  density  wake  is  considerably  thicker 
than  the  velocity  wake,  with  the  temperature  wake  falling  between  the 
two.  Note  that  these  differences  would  be  even  more  prominent  if  the 
Prandtl  number  had  not  been  taken  as  unity. 

Finally,  Figure  16  shows  the  physical  "velocity  half-thickness"  b* 
for  two  different  Mach  numbers  and  temperature  ratios.  The  b*  is  defined 
as  the  distance  from  the  wake  plane  of  symmetry  and  the  point  whre  u  is 
0.01,  and  was  computed  using  a  transformation  inverse  to  that  of  eq.  (16). 

We  see  that  (a)  the  wall  temperature  does  not  much  affect  the  low-speed  wake 
thickness,  but  it  does  so  for  =  6,  (b)  the  high-speed  wake  is  thicker 
by  a  factor  of  2  to  4  than  the  low-speed  one  depending  on  T  /T  (this  is 
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really  the  effect  of  the  initial  S /fl i)  and  (c)  the  wake  grows  by  only  a 
small  amount  (about  a  factor  of  2  at  most)  before  x'  becomes  unity. 

The  above  results  illustrate  the  findings  of  Section  5  that  the  present 
theory  achieves  two  major  goals:  an  uneventful  transition  from  the 
discontinuous  T.E.  profile,  and  agreement  with  the  classical  asymptotic 
behavior  at  large  x'.  Regarding  the  latter,  the  early  papers  (Ref.  7) 
indicate  that  the  asymptotic  solution  is  reached  within  about  "three 
plate  lengths"  (i.e.  x*  =  31)  of  the  hypothetical  flat  plate  of  length  l 
generating  the  wake.  We  can  check  to  see  if  the  present  theory  agrees  by 
connecting  our  coordinate  x'  with  the  equivalent  normalized  distance  x*/l. 
Using  the  customary  incompressible  formula  for  the  T.E.  0  in  terms  of  l: 

=  (e.twfjL 


we  get,  from  eq.  (36): 

v  ▼  yr  X 

f  y*  c  X_.  _  „  *  .  _  _ _  v  -  0,5&  — 

X  S  40Rfe  "  4&l?ee  ff.4W)fc(.  I 

At  x*/£  =  3  or  x'  =  1.68  our  present  theory  feq.  (47)  or  Figure  4  )  gives 
w  =  0.206,  while  the  asymptotic  value  from  eq.  (48)  is  0.218  for  a 
difference  of  only  about  5  percent;  thus  the  three-plate-length  estimate 
still  holds.  On  the  other  hand,  the  corresponding  comparison  fo  -he 
density  defect  is  0.345  and  0.589  for  the  actual  and  asymptotic  values 
respectively,  for  a  difference  of  70  percent  (assuming  M  =  5  and  TW/TQ  =  0.2). 

7.  Computations  and  Graphic  Results  for  the  Asymmetric  Wake 


As  already  noted,  the  computations  for  asymmetry  were  limited  to  the 

extreme  case  of  P  =  0.  To  keep  this  report  to  a  reasonable  size,  numerical 

computations  with  eqs.  (71),  (72),  etc.  were  done *for  three  cases:  Me  *  6, 

T  =  0.2T 
w  o 


Me  =  °.  Tw  s  °’2To;  and  Me  “  6*  Tw  “  V 


Comparison  of  the 

former  two  a  lows  an  appreciation  of  the  compressibility  effect,  while  the 

heat-transfjr  effect  can  be  illustrated  to  some  extent  by  comparing  the 

first  with  the  third  cases.  Of  course,  the  results  for  the  velocity  are 

valid  for  all  M  and  T  ,/T  because  of  the  transformed  coordinates, 
e  w  o 


A  general  illustration  of  the  course  of  events  beyond  the  trailing 
edge  is  given  by  Figures  14  through  21  which  include  the  initial  variations 
(at  the  T.E.)  of  the  variable  plotted.  Note  that  these  initial  variations 
are  uniformly  zero  for  all  y'  >  0,  since  P  =  0  implies  that  01  =  0  if 
@2  is  finite;  that  is,  there  is  no  boundary  layer,  initially,  on  the 
positive  half-plane.  Also  note  that  in  these  Figures,  the  plotted  quantities 
are  the  variations  i  -  u/ue  etc.  rather  than  the  profiles  as  defined  by 
eq.  (75)  etc. 


Just  as  in  the  symmetric  wake,  the  discontinuity  vanishes  as  soon  as 
the  fluid  leaves  the  T.E.,  the  "rounded  top"  of  each  variation  marking  the 
"wake  center".  The  T.E.  asymmetry  persists  for  some  distance,  symmetr'zation 
occurring  gradually  toward  the  larger  x'.  The  wake  center  deflects  into 
the  y'  <  0  half  plane;  for  the  velocity,  which  always  starts  out  with  the 
center  at  y'  =  0  the  deflection  Is  gradual.  For  the  temperature,  the  Mg  a  6, 
T  =  0.2T  case  shown  in  Figure  20,  the  center  initially  lies  off  the 
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y'  *  0  plane,  l.e.  at  Infinite  x'  the  wake  symmetrizes  but  remains  displaced 
off  the  y'  »  0  plane.  For  cases  of  practical  Interest,  as  Figure  20  shows, 
the  deflection  does  not  amount  to  more  than  one  or  two  .  However,  the 
fact  that  the  wake  center  for  the  velocity,  called  y'  ,  does  not  coincide 
with  that  for  the  temperature  or  density  (y*  ,  y'  respectively)  Is  a 

v  V 

major  finding  of  this  work.  Figure  22,  on  which  the  deflections  of  y'  , 
y'r  *  and  y'  are  drawn,  shows  that  the  centers  for  these  three  properties 
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remain  deflected  along  different  path.,  as  x'  becomes  Infinite.  This 
difference  In  the  wake  center  of  velocity,  temperature,  and  density  lies 
In  the  fact  that  whereas  the  velocity  center  is  always  at  y'  =  0  when 
x'  =  0,  the  extremum  in  the  temperature  distribution  at  the  T.E.  may  be 
off  the  y'  =  0  plane  (the  case  of  the  cooled  hypersonic  wake  discussed 
here  is  a  typical  example). 

The  above  remarks  refer  to  the  compress ible-transformed  deflection; 
in  the  physical  plane,  the  wake  centers  (for  the  velocity  only)  are  shown 
on  Figures  23,  24,  and  25.  It  was  found  that  in  the  physical  plane,  the 
deflection  decreases  at  large  x'  generally,  l.e.  the  wake  centers  tend  to 
return  to  y'  *  0.  Since  the  deflection  amounts  to  only  a  small  physical 
distance  anyway,  it  probably  deserves  little  additional  discussion.  Finally, 
in  the  non-dimensional  coordinate  F  defined  by  eq.  (50),  the  wake  centers 
for  all  variables  seem  to  always  return  to  y'  -  0. 

Figures  26,  27,  and  28  show  the  net  result  of  the  asymmetry  on  the 
maxima  or  minima  of  the  flow  variables  at  the  wake  center.  This  is  ex¬ 
pressed  in  terms  of  the  defects  such  as  of  eq.  (44),  (51),  etc.  except 
that  the  center  property,  being  off-axis,  is  designated  by  the  subscript 
"c".  Dashed  lines  on  these  Figures  represent  the  symmetric-wake  defects 
from  eq.  (44),  (51),  etc.  for  comparison.  It  is  seen  that  even  the  largest 
asymmetry  does  not  produce  changes  in  the  defects  of  any  more  than  a  few 
percent,  and  these  changes  naturally  vanish  at  large  x'. 

Figure  29  shows  examples  of  the  developing  profile,  while  Figure  30 
shows  how  fast  the  asymmetric  profiles  approach  the  asymptotic  limit 
expressed  by  eq.  (49).  The  example  chosen  Is  M  =  6,  T ■  0.2T  and 
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x  •  1,  that  Is  near  the  region  where  the  symmetric  wake  itself  (P  *  1) 

was  beginning  to  coincide  with  the  asymptotic  profile.  The  ordinate 

is  normalized  wUh  the  center  deficit  or  T  -  T  ,  while  the 

e  c  c  e 

abscissa  F  -  Fc  where  , 


is  shifted  to  put  the  top  of  the  profile  at  zero.  It  is  seen  that  at 
x"*  ■  1 ,  only  slight  evidence  of  the  asymmetry  still  remains,  and  that 
all  profiles  regardless  of  P  value  are  nearly  the  same  as  the  asymptotic 
one. 


Finally,  Figures  23,  24,  and  25  show  the  geographical  features  of  the 
Mg  =  6,  Tw  *  0.2Tq  wake  In  the  physical  plane,  consisting  of  the  wake  edges 
b*/&  in  the  two  half  planes  and  the  physical  wake  center  y*.  The  edges 
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and  the  center  were  again  computed  by  the  transformation  inverse  to  (16),  and 
are  limited  to  the  velocity  field;  that  Is  the  edge  Is  the  position 
where  the  velocity  profile  Is  0.01,  while  the  center  is  the  maximum  of 
the  velocity  profile  (one  could  equally  well  draw  analogous  curves  for 
the  temperature  and  the  density).  The  following  are  obvious  from  these 
Figures:  (a)  the  initially  zero-thickness  portion  of  the  wake  grows 
rapidly,  while  the  initially  finite  half  (the  one  with  initial  momentum 
thickness  02)  grows  slowly  and  on  occasion  decreases  slightly  In  thick¬ 
ness  at  the  beginning,  (b)  the  wake  center  originally  deflects  downwards 
and  then  returns  toward  the  plane  y'  s  0,  (c)  the  asymmetry  decreases 
downstream  since  the  two  edges  begin  becoming  equidistant  and  (d)  the 
wake  thickness  in  terms  of©  is  independent  of  P  at. large  x'. 

8.  Discussion  and  Comparison  With  Experiments 

The  present  work  was  motivated  largely  by  the  desire  to  provide 
closed-form  solutions  for  the  symmetric  wake  for  application  to  gas-dynamic 
lasers.  In  an  earlier  paper  (Reference  8  )  the  present  author  formulated 
a  theoretical  compendium  of  wake  predictions  for  compressible  two-dimensional 
wakes  at  arbrltrary  and  TW/TQ.  This  earlier  theory  was  based  on  an 
Integral  approximation  to  the  total  drag  of  and  heat  loss  from  the  structure 
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separating  the  two  streams,  as  well  as  the  asymptotic  symmetric  wake 
formulas  contained  herein  (eq.  (48),  (49),  etc.)  Experiments  by  Peterson 
(Reference  9  )  and  this  author  (Reference  10)  produced  data  agreeing 
with  the  turbulent-wake  prediction  of  this  earlier  theory  but  disagreeing 
with  Its  laminar  counterpart.  This  could  be  explained  by  the  fact  the 
turbulent  wake  mixes  so  quickly  that  Its  non-equl librium  (adjustment) 
zone  behind  the  T.E.  never  plays  a  role;  for  the  laminar  wake,  however, 
the  non-equilibrium  zone  (found  herein  to  extend  as  far  as  x'  «  1  to 
*  3)  was  unaccounted  in  the  early  theory,  thus  producing  disagreement 
with  the  data.  It  therefore  became  Important  to  have  on  hand  a  wake 
theory  treating  equally  the  asymptotic,  far  region  as  well  as  the  non¬ 
equilibrium  ("near"  or  "adjustment")  region  close  to  the  T.E. 

The  same  data  of  Reference  10  disagreeing  with  the  earlier  theory 
of  Reference  8  are  plotted  on  Figure  5  .  Two  data  points  from  Sato 
and  Kuriki  (Reference  11)  are  also  included.  The  agreement  Is  very  good 
within  the  data  scatter.  No  other  quantitative  data  could  be  found  for 
comparison,  but  neither  could  intimations  be  found  In  the  literature 
contradicting  the  present  findings.  In  fact,  Batt  and  Kubota  (Reference 
12)  present  an  experiment  frequently  intimating  the  validity  of  the  present 
approach;  for  example,  they  find  that  the  proper  scaling  length  is  the 
total  T.E.  momentum  thickness. 

Both  Figures  4  and  5  show  how,  in  the  limit  Me  =  0,  Tw  s  TQ,  the 
present  theory  agrees  with  the  findings  of  Toll  mien  and  Goldstein  for  the 
entire  x"  range.  This  agreement,  as  well  as  the  agreement  with  the  data, 
refers  to  the  symmetric  wake.  No  previous  theoretical  or  experimental 
threatment  of  the  asymmetric  wake  could  be  found. 

Note  should  be  made  of  the  differences  between  the  adjustment  of  the 
density  field,  and  the  adjustment  of  the  other  flow  variables.  The  density 
is  generally  very  slow  to  adjust  to  the  asymptotic  limit.  The  reader  can 
expand  on  this  statement  by  making  numerical  computations  with  the  programs 
listed  in  the  Appendix. 
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The  good  analytic  behavior  of  the  solutions  presented  here,  and  the 
agreement  with  the  data  appear,  at  first  glance,  to  be  disproportionate 
to  the  severity  of  the  basic  assumptions  (10)  and  (11).  Recall,  however, 
that  the  solutions  are  forced  to  obey  the  correct  boundary  (Initial) 
conditions  at  the  T.E.,  while  the  method  supplied  by  Gold  guarantees  the 
proper  asymptotic  behavior  as  well.  With  the  solutions  thus  "tied 
down  to  the  proper  limits  at  the  two  ends",  It  Is  not  surprising  that 
the  Intermediate  regions  give  the  proper  results. 

9.  Conclusions 


1.  Closed-form,  analytic  solutions  have  been  found  for  steady, 
laminar  two-dimensional  wakes  applicable  for  any  value  of  the 
edge  Mach  number  M^,  tral ling-edge,  wall-to-stagnatlon  temperature 
ratio  Tw/Tq  and  for  any  ratio  P  of  the  momentum  thickness  of 

the  two  merging  boundary  layers.  The  streamwise  development  of 
the  wake  is  controlled  by  the  total  momentum  thickness  and  the 
unit  Reynolds  number  contained  In  the  longitudinal  variable  x'. 

2.  Beyond  a  certain  distance  x'  ranging  from  1  to  3  the  wake  conforms 
to  the  classic  similarity  form  regardless  of  Mg,  TW/TQ  or  the 
Initial  asymmetry  expressed  by  P.  This  distance  marking  the 
onset  of  dynamic  "equilibrium"  1?;  almost  exactly  equal  to  that 
found  in  early  studies  with  Incompressible  flows.  Only  the 
density  lags  behind  the  other  flow  variables  In  reaching  equili¬ 
brium  within  this  distance. 

3.  Many,  if  not  most,  current  practical  applications  involve  wakes 
within  the  "non-equilibrium"  distance  between  the  trailing  edge 
and  x'  =  1  to  3.  The  classic  asymptotic  formulas  are  Inapplicable 
in  this  region.  For  wakes  such  as  encountered  In  laser  cavities 
(large  M  ,  T  <  <  T  )  the  present  theory  reveals  many  peculiarities 
of  behavior,  such  as  non-monotonic  development,  off-axis  maxima 
etc.  even  for  the  initially  symmetric  wake. 
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Even  for  the  largest  Initial  asymmetry,  the  flow  defects  (e.g. 
minimum  velocity)  are  almost  equal  to  those  of  the  symmetric 
wake.  In  the  non-equlllbrljijw-'reglon,  however,  the  Initial 
asymmetry  causes  the. ..watt'" to  deflect  towards  the  side  of  the 
initially  thicker  boundary  layer,  and  the  extrema  in  the  velocity 
profile  do  not  coincide  with  the  extrema  for  the  density  or 
temperature. 
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Figure  7:  Temperature  decrement  for 
cooled  symmetric  wakes. 
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Figure  9:  Temperature  profiles  for 
the  low-speed,  cooled 
symmetric  wake. 


Figure  10:  Temperature  profiles  for 
the  hypersonic,  adiabatic 
symmetric  wake. 
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Figure  12:  Density  defects  for  the 
symmetric  wake  (typical). 
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wake. 


defects  for  the  supersonic 
cooled  symmetric  wake. 


Figure  17:  Development  of  the  velocitj 
variation  for  the  asym¬ 
metric  wake  past  the  T.E. 


Figure  18:  Temperature  development  of 
the  low-speed,  cooled 
asymmetric  wake  past  the 


Figure  19:  Temperature  development  of 
the  hypersonic  adiabatic 
asymmetric  wake  past  the 


igure  20:  Temperature  development  of 
the  hypersonic  cooled 
asynmetric  wake  past  the 


trie  wake  past  the 
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Figure  23:  "Map"  of  the  low-speed 
cooled  asymmetric  wake 


igure  24:  "Map"  of  the  hypersonic 
adiabatic  synmetric  wake 
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Figure  27:  Effect  of  asymmetry  on 
the  temperature  defect. 


Figure  29:  Comparison  of  normalized 
profiles  for  the  hyper¬ 
sonic,  cooled  asymmetric 
wake  in  the  non-equili¬ 
brium  zone. 


Figure  30:  Decay  of  initial  profile 
asymmetry  In  the  hyper¬ 
sonic,  cooled  wake. 
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APPENDIX 

COMPUTER  PROGRAMS 


All  programs  contained  in  this  Appendix  use  the  BASIC  language. 

The  output  symbols  for  each  program*  as  they  appear  on  the  printout* 
are  explained  in  terms  of  the  symbols  listed  in  the  beginning  of  this 
report. 


> 
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WAKEFL0W2 "  PROGRAM 


Purpose:  To  compute  the  defects  for  symmetric  wakes  as  a  function  of 
x",  with  Mg  and  TW/TQ  as  parameters. 

Application  Example 

•  '  i  • 

2.1  If  the  defects  at  a  single  x',  Mg,  TW/TQ  are  desired: 

Compute  the  defects  at  M  =  3,  T  /T  =  1.5,  x'  =  0.6. 

c  WO 

Enter: 

30  LET  T  =  1.5 
40  ,  LET  M  -  3 
180  LET  x  -  0.6 
410  -* 

460 

470 

2.?  If  tho  defects  at  a  single  Mg,  T(^/T0  but  a  series  of  x'  are  desired 
Compute  the  defects  for  the  above  example  but  for  x"  intervals 
o*  0.05  from  0  to  0.8. 

Enter: 

30  LET  T  =  1.5 

40  LET  M  =  3 

180  FOR  x  =  0  TO  0.8  STEP  0.05 

460 

470 

2.3  If  one  needs  the  defects  for  a  series  of  x*  at  each  of  a  series 
of  Mg,  and  additionally  for  a  series  of  Tw/To: 

Compute  the  defects  for  M  =  0,  1 ,  2.  .  .  5,  and  T  /T  =0.5, 

1,  1.5,  each  for  x'  from  0  to  2  in  steps  of  0.2. 

Enter: 

30  FOR  T  =  0.5  TO  1.5  STEP  0.5 
40  FOR  M  »  0  TO  5  STEP  1 

180  FOR  x  =  0  TO  2  STEP  0.1 
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3.  Explanation  of  Output  Symbols 


Printout  Symbol  Explanation 


B  B 

C  C 

TO/TE-1  t 

TO-TE/TW-TE  V 

1-R0/RE  r 

RE-RO/RE-RW  r 

W  w 

M{0)  M(0) 

T00-T0E/T0E  H 

TOO-TOE/TW-TOE  H 
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"WAKEFL02"  LISTING 


REM  WAKE  FI.  02  !  . 

DXff  X<40  )  , X  I.  <40  >  ,Y<40>  »W<40) ,X2<40  ),HHO  )  ,f'MO)  ,2(4/0 
DX'fl  H1(407,R<40>»R1<40>  ■ 

DIM  ,ri;t  (40  )  j.  H2  ( 4 0  >  ,  H3 <  *i {) ) 


FOR  T*. 
FOR  M 
PRINT  1 
PR  IN  I'  ’ 
PRINT  1 
PRINT 
J-  LEI  B= 
1  LET  CL 
]:.  PRINT 


,  2'  TO  .1,0  STEP  .  4 
0  TO  5  STEP  1 
"SYMMETRIC  WAKE  DEFECTS'" 
"FOR  TW/TO'="M 
"AND  MACH  NO-'"  tM 


L.E1  B” » 2*  (MA2 )  *  ( T+l )  +  T--1 
LET  C-  -.2*(hA2> 

PRINT  "THE  CONSTANT  B*"*© 

PRINT  "THE  CONSTANT  C>  "  J  C 
PRINT 

PRINT  "TO/TE-l  "  ,  "TO  -TE/TW-TE"  ,  "1 -RO/EE"  ,  "EE--RO/T, 

PRINT  TAB(R)  ,  "W"  , "MOO  > ",  "TOG-- TOE/TOE"  »  "TO  O-TOE/TW-TOE 


E--RW" 


5  PRINT  TAP  CO) , "W" , "M( 0 ) " ,  "TOO-  TOE/TOE" , "TO  0 -T OE/TW-TOE 

y,  PRINT  ■- 

)  FOR  X  “ 0  TO  3  STEP  .3 

/  LET  Pit* 327591 

)'  LET  Al i::: .  7.  54329092  il 

)  LET  A  2  ?B4 19673600 0  0 00  j. 4 

)  LET  A3«t. 421413741* 

)  LET  A4«~l, 453152027# 

)  LET  A5-1 . 061405429# 

)  LET  X1«2#SQR<X> 

)  LET  D«l/.<  l't-p)KXl) 

)  LET  Y«< A1*D+A2*<DA2>+A3*<DA3> fA4*<DA4)+A5*<DA5> >#EXP( 

1  LET  W«<\EXP<4#X)  )*Y  . 

!  LET  .X2L2*X1 
)  LET  F--:l/<  1  t-P*X2) 

)  L  E  T  2  “  ( A 1  *  F + A  2  *  <  F  A  2 )  +  A  3  *  ( F  A  3 )  +  A  4  *  ( F  A  4  )  +  A  5  *  <  f  • ' A  5 )  >  *  E  X  P  <  ■ 
)  L E T  H - B * W  +  C * E X P  <  1. 6*X>  * 7. 

)  LET  m-»M/(B+C) 

1  LET  R'-  H/  (H+l ) 

)  LET  R:L"R*  < :L  +  1/  <  B •♦■(:: ) ) 

)  LET  Ml®h*(l--W)/SQR<1+H) 

)  L  E  T  H  2  -  (  ( H  + 1 )  ( J.  .  2  *  ( M  .1 A  2 )  )  /  ( 1  .  2  *  ( M A  2 )  )  )  -•  1 
1  LET  H3«H2/<T-1? 

)  PRINT  X »  H , H 1 1 R  f E 1 
)  P  R I N  T  T  A  F:  (  8  >  ,  W  ,  M 1 ,  H  2  ,  H  3 
1  NEXT  X 
)  PRINT 
)  PRINT 
)  NEXT  H 
)  NEXT  T 
)  END 


(XI A2 )  ) 


■  ( X2A2 )  ) 


LET  H2= 
LET  H3: 


«^***:~  ^4)**;* , -f --  -fill  ''*•»**• 


"WAKEFL02"  OUTPUT 


3  Y  M  M  F  T  R I  c  W  A  K  E  D  E  F  E  C  T  S 

FUR 

TW/TO  --  .5 

AND 

MACH  NO-  2 

THE 

CONSTANT  B-  .7 

THE 

CONSTANT  C- ••♦8 

V  / 

TO/TE-1 

TQ-TE/TW--TE 

1-R 0 /RE 

RE-RO/RE-RW 

W 

M<  0  ) 

TOO -TO E/TOE 

TOO-TOE/TM-TOE 

0 

•••• .  :l. 

1. 

-  ♦  1  Jill 1 

1 

1 

0 

•“ » 5 

1 

.2 

*  0967988 

•••• ,  967988 

♦  0382557 

-.794302 

♦438246 

.1.  ♦  0  3459 

-.260224 

♦520449 

.  4 

*  0883864 

- . 883865 

.0812087 

- ♦ 730879 

♦  364732 

1 ♦ 2 1 785 

-.215979 

. 431957 

*  6 

.  080  449.1. 

••■•♦80  4491 

♦074459 

•v  670 131 

♦  3.1.37.I6 

1.32048 

-.190423 

.380847 

.  8 

.  0740363 

••••  .74  0  563 

.0689501 

-♦62.9551 

♦  279967 

.1. .  38953 

•♦172881 

.345762  • 

.1 

♦  0688763 

- . 688763 

♦  064438 

-.579943 

.255403 

1 .44042 

-♦159769 

.319587 

1.  .  2 

♦  0645929 

•  ♦ 645929 

♦0606738 

•••♦546064 

♦236448 

1.48005 

-.149444 

♦298887 

:l.  i  4 

.0609823 

■•♦609823 

.0574772 

-.517295 

♦  22:1232 

1.51211 

-.141019 

.282037 

1.6 

.  0578878 

"".578878 

.0547201 

-.492481 

♦208658 

1.53877 

-.133963 

.267926 

.1.  *  8 

♦  055.1  987 

••■♦551988 

.0523112 

-.47030.1. 

. 198037 

1.56141 

-.127937 

.255873 
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"WAKEFL03"  PROGRAM 


Purpose:  To  compute  the  lateral  profiles  of  symmetric  wakes  at  a 
fixed  M  ,  T  /T  ,  a6  a  function  of  y'. 

C  WO 

Application  Example 

Compute  the  profiles  for  x'  =  0.25,  M£  =  1  and  TW/TQ  «  1,  for  y' 
intervals  of  0.15  from  0  to  3. 

Enter: 

30  LET  M  =  1 

40  LET  T  =  1 

50  LET  x  =  0.25 

3f>0  FOR  V  =  0  TO  3  STEP  0.15 

Explanation  of  Output  Symbols 

Printout  Symbol  Explanation 


Y'  y' 
i't 2SQRX'  F 
UTILDA  u 
TTILDA  T 
RH0TILDA 

M  M 
TOTILDA  H 
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"WAKEFL03"  LISTING 


1.  0 

REM  WAK 

:i.  5 

DIM  Ml  ( 

2  0 

DIM  V<2 

2 1 

DIM  R5< 

30 

LET  M-l 

4  0 

...FT  T  - 1 

5  0 

LET  X-~ . 

.f»0 

PRINT  " 

70 

PRINT  " 

80 

"RINT 

9  0 

PRINT  " 

:i.  0  o 

PRINT 

:l  1 0 

PRINT 

120 

PRINT 

130 

LET  B  •' 

140 

LET  C= 

.1.50 

PRINT 

160 

PRINT 

170 

LET  P::: 

1 80 

LET  A.l 

190 

LET  A 2 

20  0 

LET  A3 

210 

LET  A 4 

220 

LET  A 5 

230 

LET  XI 

240 

LET  0“ 

250 

LET  Y» 

260 

LET  W» 

270 

PRINT 

280 

LET  X2 

290 

LET  F  = 

30  0 

LET  Z« 

310 

LET  H* 

320 

PRINT 

325 

PRINT 

330 

PRINT 

34  0 

PRINT 

345 

PRINT 

350 

PRINT 

360 

FOR  V* 

37  0 

LET  Y 1 

380 

LET  Y2 

390 

LET  R= 

40  0 

LET  Y3 

405 

IF  Y3< 

410 

LET  Y4 

420 

LET  R'l 

421 

GOTO  4 

422 

LET  Y3 

424 

LET  Y 4 

I.  03 


FOR  X'; 
"AND  M“ 


"  JX 

"  ;m 


;e: 

JO 


"AND  TW/TO:~"  }  T 

♦  2* ( MA2 ) *  <  T+l ) +T-1 
” . 2* ( MA2 ) 

"COMPUTED  CONSTANT  B 
"COMPUTED  CONSTANT  C 

♦  327591 
*•254829572* 

—  ♦  2844967360 0  0 0  001# 

=1.421413741# 

s-1 .453152027# 

-- 1 . 061405429# 


<  A1  *D+A2*  ( DA2  )  +  A3*  (  DA3 )  +A4*  ( DA4 )  •♦•A3*  <  DAS )  )  *EXP  <  --  (  XI A2 )  ) 

Y*EXP <  4*X ) 

"COMPUTED  W""JW 

:-'2*Xl 

=1/<1+P*X2) 

» < A1 *F  +  A2* ( FA2 >  + A3* <  FA3 ) +A4# ( FA4 ) +A5* ( FAS ) ) *EXP < - ( X2A2 ) ) 
»B*W+C*EXP< 16*X>*Z 
"COMPUTED  H-T<0)/TE  -  1*"JH 

"COMPUTED  DEFECT  H '  =  ( T  <  0  )  -  TE/TN--TE ) «"  J 1-1/  ( B+C ) 

"Y '  "  ,  "Y '  /2SORX 7  "  ,  "UTILDA"  ,  "TTII...DA"  ,  "RHOTILDA" 

TAB  (  8 )  t  "M"  ,  "T0TTI...DA" 

0  TO  3  STEP  .15 
»V/X1+X1 
~ 1 / < 1+PKY1 ) 

( A1 XY2+A2* ( Y  2  A  2 ) +A3* < Y2A3 >  +  A4*  <  Y2A4 >  +A5* <  Y2A5 ) ) *EXP < - <  Y 1 A2 ) ) 
-“D/Xl-X  J. 

0  GOTO  422 
“1/ < 1+P*Y3  > 

=  1  -  ( A1XY4+A2* <  Y4A2 ) +A3* ( Y4A3  >  +A4# ( Y4 M ) +A5* ( Y4A5 ) ) *EXP < - ( Y3A2 ) ) 
30 

“Xl-V/Xl 
=  l/<l-t-P*Y3) 


« 
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"WAKEFL03"  LISTING  (CONT'D) 


126  LET  Rl~:  <  A1#Y1+A2*  <  Y1A2  )+A3* <  Y1A3 ) +A1* ( Y1A1 H  A5* ( Y1A5 ) ) #EXP  < - <  Y3A2 )  ) 

128  LET  U«.5*<R*EXP(4*X+2*V>+Ri*EXPM*X-2*V>>/N 

129  GOTO  110 

130  LET  U«.5*(R*EXP<1#X+2*V>  +  <l4-Rl>*EXP(1*X-2*V>  )/W 
110  LET  Y5®2#X1+V/X1 

150  LET  Y6~l/( 1+P*Y5) 

160  LET  R2“ ( A lw YA+A2* < Y6A2 ) +A3* ( Y6A3 ) +A1* ( Y6A1 ) +A5* < YAA5 ) ) #EXP ( - ( Y5A2 ) ) 
170  LET  Y7--X1-V/X1 
180  IF  Y7<0  GOTO  520 
190  LET  Y8=1/<1+P*Y7> 

500  LET  R3« ( A1*Y8+A2* <Y8A2)+A3*(Y8A3)+A1*  Y8A1 ) +A5* < Y8A5 > > KEXP ( - ( Y7A2 ) > 

510  GOTO  550 

520  LET  Y7-0/X1 -Xl 

530  LET  Y8- 1/ < 1+P*Y7 ) 

510  LET  R3=2~<A1*Y8+A2#<Y8A2)+A3K<Y8A3)+A1*<Y8A1)+A5*<Y8A5) ) *EXP ( - < Y7A2 ) ) 
550  LET  Y9»2*X1-V/X1 
560  IF  Y9<0  GOTO  600 
57  0  L ET T  Fl-"l/<  l*fF:,#Y9> 

500  LET  R1=  ( A1*F1+A2# <  F  1 A2 )  +A3*(  FI  A3 )  +A1*(  E.1  A1 )  +A5*(F.l A5 ) ) *EXP ( ~ ( Y9A2 )  ) 

590  GOTO  630 

600  LET  Y9«V/Xi-2*Xl 

610  LET  FI  *l/<j>P*Y9> 

620  LET  R1*2~<Al#Fl+A2*<FlA2>+A3«('FiA3>+A1#<FlA1)+A5*<FlA5> >*EXP<-< Y9a2> ) 
630  LET  Tl * ♦  5* ( E!* (EXP <1*X-2*V ) *R3+EXP <4*X+2*V >*R>> 

633  LET  T2“  •  5*  < C#  < EXP ( 16*X~1*V )  *R1+EXP ( 1 6*X+1*V )  *R2  >  ) 

637  LET  T 3 ••••TUT 2 
610  LET  T 1--T3/H 

615  LET  R5«  (  Tl*  <  H+.1 )  )  /  <  T1*H+1 ) 

616  LET  M1~M#  <  1--U*W )  /SQR  < 1+T1#H) 

617  LET  H2» ( < H+l ) * ( 1  + . 2* ( hi A2 ) > / ( 1 ♦ .  2*  <  MA2  > ) ) -1 

618  LET  H3“ ( 1  +T4*H )  #  <1  + .  2*  <  Ml  A2  >  >  /  <  1  +  ♦  2*  < MA2 )  > --1 

619  LET  H1-H3/H2 

650  PRINT  V  t D/Xl t U >  T1»R5 
655  PRINT  TAB ( 8 ) , Ml , HI 
660  NEXT  U 
700  END 
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"WAKEFL03"  OUTPUT 


SYMMETRIC  WAKE  PROFILES 


FOR  X ' ~  .25 
AND  M~  1 
AND  TW/T0«:  1 

COMPUTED  CONSTANT  B»  .4 

COMPUTED  CONSTANT  C«-.2 

COMPUTED  W:s:  .427 584 

COMPUTED  H~T ( 0 ) /TE  -  i»  .119953 

COMPUTED  DEFECT  H' - <  T <  0 ) -TE/TW-TE ) =  .599765 


Y  ' 

Y ' /2SQRX ' 

M 

UTIL DA 

T0TILDA 

TTILDA 

RHOTILDA 

0 

0 

.540894 

1 

1 

1 

1 

.  15 

.15 

.54702 

.985738 

1  .  12106 

.987599 

.988913 

,3 

.3 

.564875 

.944288 

1.7231 

♦951407 

.956384 

.  45 

.45 

. 592966 

.879449 

11.5166 

.894323 

.904561 

.6 

♦  6> 

♦629039 

.796859 

-1.70369 

.820713 

♦836781 

.75 

.75 

.670389 

♦703141 
- .675035 

.73582 

.757247 

♦  9 

.9 

.714203 

♦604977 

-.371414 

.645114 

.670604 

1.05 

1.05 

.757806 

.508323 

-.227826 

.553693 

.58149 

1 .2 

1.2 

.799303 

♦417857 

-.145981 

.465832 

.4941 

1.35 

1.35 

.836925 

.336733 

-.0949776 

.384717 

.411859 

1.5 

1.5 

.869852 

.266615 

-.0618271 

.312385 

.337221 

1,65 

1.65 

.897756 

.207897 

-.03993 

.249804 

.271629 

1.8 

1.8 

♦920752 

.160037 

-.0254572 

. 197067 

.215609 

1.95 

1.95 

.93926 

.  121906 
-.0159722 

.153633 

.168949 

2.1 

2.1 

.953864 

.0920926 

-9.84869E-03 

.118563 

.130922 

2.25 

2.25 

.965202 

♦0691332 

-S.96667E-03 

.0907185 

.100507 

2,4 

2.4 

♦973094 

.0516595 

-3.55195E-03 

.0689246 

.0765593 

4.1  ♦  \J \J 

2 . 55 
.98049 

.0384782 
-2.08039E -03 

,0520665 

.05795 

2 . 7 

2.7 

.98546 

. 0285984 
-1.20079E-03 

.  0391517 

.0436431 

2,85 

2 . 85 
.989185 

.0212262 
-6.84796E- 04 

. 0293344 

.032738 

65 


"WAKEFL04"  PROGRAM 


1.  Purpose :  To  compute  the  physical  velocity  half- thickness  (defined  at 

u  »  0.01)  for  a  symmetric  wake,  given  Me,  Tw/TQ,  x',  and 
y'(u  *  0.01).  The  latter  Is  found  by  first  using  the 
WAKEFL03  program. 

2.  Application  Example 

Compute  the  half  thickness  for  Mfi  *  5,  yTQ  =  1,  x'  *  0.1;  for  this 
case,  WAKEFL03  gives  y'(u  *  0.01)  *  2.7. 

Enter: 

30  LET  M  =  5 

40  LET  T  *  1 

50  LET  X  =  0.1 

360  LET  V  *  2.7 
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"WAKEFL04"  LISTING 


10  REM  WAKEFL04 

2  0  DIM  U(40>  ,Y1.  (40)  ,Y2<40>  »R<40>  ,Y3<40>  , Y4< 40 > , R1 < 40 )  ,U<40> 

21  DIM  R5<4U),V1(40>,R6<40>,V2<40> 

3  0  LET  M*5 

4  0  LET  T*1 
50  LET  X».l 

60  PRINT  "PHYSICAL  VELOCITY  HALE-THICKNESS  COMPUTATION  (LJTILDA -0 . 01 )  " 

70  PRINT  "FOR  SYMMETRIC  WAKES" 

72  PRINT  " . — - - - " 

SO  PRINT 

90  PRINT  "FOR  X'«"$X 
1.0  0  PRINT  "AND  M*"JM 
110  PRINT  "AND  TH/T0«"?T 
120  PRINT 

130  LET  B- *  2* ( MA2 ) *  <  T+l )  +  T~1 
140  LET  C~~ ♦ 2* < MA2 ) 

170  LET  P-. 327591 
180  LET  Al»*254829592# 

190  LET  A2*~. 2844967360000001* 

200  LET  A3»l ♦421413741# 

210  LET  A4*-l ♦453152027# 

220  LET  A5-1 • 061405429* 

230  LET  X1»2*SQR<X> 

240  LET  D«1/<1+P*X1> 

25 0  LET  Y* <  A 1 *D+ A2* < DA2 ) + A3* ( DA3 ) +A4* ( DA4 ) + AS* ( DA5 ) ) *EXP ( - <  XI A2 ) ) 

260  LET  W«Y*EXP<4*X> 

28 0  LET  X2“2#X1 
290  LET  F-l/  <  1+P*X2  ) 

30  0  LET  Z- ( A1*F+A2* ( FA2 ) +A3* ( FA3 ) +A4*< FA4 ) +A5* < FA5 ) ) *EXP < ~ ( X2A2 ) ) 

310  LET  H-B*W+C*EXP<1A*X>*Z 

360  LET  V«2.7 

650  FOR  I«1  TO  40 

660  LET  VI <I)“I*V/40 

665  LET  Vl-VKI) 

1370  LET  Y1-V1/X1+X1 
1380  LET  Y2»1/<1+P*Y1> 

1390  LET  RMA.t#Y2+A2#<Y2A2)+A3#<Y2A3> +A4* < Y2A4 ) +A5# <Y2A5) )*EXP<-( Y1A2) ) 
1400  LET  Y3-V1/X1-X1 
1405  IF  Y3<0  GOTO  1422 
1410  LET  Y4*1/(1+P*Y3) 

1420  LET  Rls 1-<A1*Y4+A2*< Y4A2 ) +A3# < Y4A3 > +A4# < Y4A4 > +A5# < Y4A5 > > *EXP< - ( Y3A2 ) > 

1421  GOTO  1430 

1422  LET  Y3-X1-V1/X1 
1424  LET  Y4:;::  1/(1  +P*Y3 ) 

1426  LET  R1*(A1JKY4+A2*<Y4A2)+A3*<  Y4A3>+A4*<  Y4A4)+A5*<  Y4A5) )  *EXP  <  -  ( Y3A2 )  ) 

1428  LET  U”  ♦  5*  (  R*EXP  (  4*X+2#V.l )  +R1*EXP<  4*X  ~2*V1 )  )  /W 

1429  GOTO  1440 

1430  LET  U- .  5*  <  R*EXP  ( 4*X+2*V1  >  +  ( 1 +R1 )  *EXP  ( 4*X~2*V1 )  )  /W 
1440  LET  Y5-2*X 1+V1/X 1 
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"WAKEFL04"  LISTING  (CONT'D) 


1450  uirr-Y6"l/  ( 1+PKY5) 

1460  LET  R2:“(A1#Y6+A2*<Y6A2MA3*<Y6A3>+A4*(Y6A4>+A5#<Y6A5) >*EXP<-< Y5A2> ) 
.1470  LET  Y7*X1“V1/X1 
1480  IF  Y7<  0  GOTO  1520 
1490  LET  Y8»1/<1+P*Y7) 

1 50 0  LET  R3» ( A1 * Y8+A2* ( Y8A2 ) +A3* < Y8A3 ) +A4*< Y8A4 > +A5# ( Y8a5 > ) *EXP ( - ( Y7A2 ) ) 

1510  GOTO  1550 

1520  LET  Y7«Vi/Xl-Xl 

1530  LET  Y8»1/<1+P*Y7) 

1540  LET  R3»2-(A1KY8+A2x< YSA2 > +A3* < Y8A3 > +A4* < Y8A4MA5*( Y8A5> >*EXP(-(Y7A2> ) 
1550  LET  Y9»2*X1-V1/X1 
1560  IF  Y9<0  GOTO  1600 
1570  LET  Fl»l/a+P*Y9) 

1 58 0  I..E T  F\' 4  -  ( A 1  *F  1 +A2*  <  F  I  A2 )  + A3*  ( F 1 A3  >  +A4*  <  FI  A4 )  +A5*  ( F 1 A5 )  )  *EXP  ( -  <  Y9 A2  > ) 

1590  GOTO  1630 

1600  LET  Y9'"0.l/Xl-2*X1 

1610  LET  F1«1/<1+P*Y9) 

1620  LET  R4*2- -< A1*F1+A2*<F1A2>+A3*<F1A3>+A4*< F1A4 >+A5*<FlA5 ) >  *EXP  <  -  <  Y9A2  ,> > 
1630  LET  T.L « ♦  5*  < B*  < EXP  ( 4*X-2*V.l )  XR3+EXP  ( 4*X+2*V1 )  *R  >  > 

1633  LET  T2« ♦  5*  <  C*  < EXP  < 16*X- A*V1 ) *R4+EXP ( 16*X+4#V1 ) *R2 ) ) 

1637  LET  T3»T1+T2 

.164  0  LET  T4=“ T3/H 

.1. 6  4  5  L.  E  T  R  5  « <  T  4  *  ( H  + 1  >  )  /  ( T  4  *  H  +  .1 ) 

1650  LET  RA«1+T3 
1 660  L  E  T  02  «  V  2 + R  6  *  O  /  4  0 
1670  NEXT  I 

180  0  PRINT  "TRANSFORMED  HALF  THICKNESS  <  YTILDA/CAPTHETA )  I V 
1810  PRINT  "PHYSICAL  HALF  THICKNESS  ( Y/CAPTHETA ) »" J 02 
2000  END 
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"WAKEFL04"  OUTPUT 


PHYSICAL  VELOCITY  HALF-THICKNESS  COMPUTATION  <UTILDA“0 ♦ 01 ) 
FOR  SYMMETRIC  WAKES 


FOR  X'“  .1 
AND  M^  5 
AND  TW/TO*  1 

TRANSFORMED  HALF  THICKNESS  ( YTILDA/CAPTHETA > 58  2.7 
PHYSICAL  HALF  THICKNESS  < Y/CAPTHETA ) *  6.29334 
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"WAKEFL06"  PROGRAM 


t 


I 

I 


1.  Purpose:  To  compute  the  "profiles"  for  the  completely  asymmetric 

(P  «  0)  wake  for  a  given  x',  Me,  and  TW/TQ,  for  intervals 
l  of  y'  over  a  specified  range  of  yA. 

2.  Application  Example 

Compute  the  profiles  for  P  *  0  and  x*»  0.01,  ■  1  and  TW/TQ  a  1.01 

l  in  steps  of  y'  *  0.3  over  the  range  from  y'  ■  -3  to  +3. 

Enter: 

"  20  LET  X  *  0.01 

22  LET  M  *  1 
24  LET  T  »  1.01 

[  160  FOR  Y  =  -3  to  3  STEP  0.3 

I 

3.  Explanation  of  Output  Symbols 

:  Printout  Symbol 


Y72SQR(X') 
1  -  U/UE 
T/TE  -  1 
R/RE  -  1 
M 

TO/TOE  -  1 


C 
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Explanation 


F 

1  _ 

VTe 

M 


u/ue 
-  1 
-  1 


VToe 


-  1 


( 


% 


"WAKEFL06"  LISTING 

.10  REM  HAKEFL06 

.15  DIM  Y<40>»Z(40)»E:<40>»Y1<40>»U<40> 

16  DIM  21  <40>  ,D<40> ,  Y2<40>  ,Y3<40)  ,H<40>  ,R<40> , M < 40 > » T1 ( 40 > 

20  LET  X«=.01 
22  LET  M™1 
24  LET  I  “ .1*01 

26  LET  ESI™  ♦£#  <MA2 )* < T+i  )  +T--1 
28  LET  Ca-.2*<MA2> 

30  LET  P». 327591 
40  LET  Al™ « 254829592# 

50  LET  A2*-* 2844967360 000001# 

AO  LET  A3»l .421413741* 

70  LET  A4™~1 . •053152027# 

80  LET  A5»l. 061405429* 

90  PRINT  "COMPLETELY  ASYMMETRIC  <P»0>  WAKE  PROFILES" 

100  PRINT  " - - - - - " 

110  PRINT 

120  PRINT  "FOR  X'“"$X 
122  PRINT  "AND  MACH  NQ.«"$M 
124  PRINT  "TW/TO*" I T 
.12 A  PRINT 

130  PRINT  "Y ' " , "Y V28QR(X# ) " ♦ "l-U/UE" » "T/TE-1" , "R/RE-1 " 

140  PRINT  TAP (8) » "M" » "TO/TQE~i " 

150  PRINT 

1 AO  FOR  Y»-3  TO  3  STEP  .3 
170  LET  Z=SQR(X>+Y/<2*SOR<X> > 

175  IF  Z<0  GOTO  192 
180  LET  B=1/(1+P*Z) 

190  LET  Yl~  <  A1*ES+A2*  (BA2  >  +A3*  <EsA3>  +A4#<BA4  )  +A5* <BA5 )  > *EXP(- ( ZA2 >  > 

191  GOTO  200 

192  LET  Z™-Y/ < 2#SQR < X > )-SGR(X) 

194  LET  B=l/<  :l+P*Z) 

19A  LET  Yi=2™< A1#B+A2*<BA2  )+A3*<BA3)+A4*<BA4)+A5#(BA5)  )*EXP<-(ZA2>  > 

200  LET  U“ ♦ 5#EXP<  X+Y ) *Y1 

210  LET  Z1~2#SQR ( X ) +Y/ ( 2*SQR ( X  > ) 

220  IF  Z1<0  GOTO  2AO 
230  LET  D=l/(  1+pxZl > 

240  LET  Y2- <  A1*D+A2# <DA2)+A3*<DA3)+A4*(DA4> +A5#  <DA5 > ) *EXP < -< Z1 A2 ) ) 
250  GOTO  290 

2A0  LET  Z1*-Y/(2*SQR<X)  )™2*SQR(X) 

270  LET  D*l/<l+PxZl> 

280  LET  Y2~2-<A1*D+A2#(DA2)+A3#<DA3)+A4#<DA4) +A5* (Da5 ))*EXP(-(ZiA2>> 
290  LET  H™  .  5*C#EXP<  4#X  +  2*Y  )  *Y2+B1*IJ 
300  LET  Y3«Y/(2*SQR<X) ) 

310  LET  R«-H/(H+1> 

320  LET  Ml™< l~U)*M*SOR< 1+R) 

330  LET  T 1  —  <  < H+l ) *<  H ♦ 2*  < M1A2 ) ) / ( 1  +  .2* ( MA2 ) ) ) -1 

350  PRINT  Y.Y3,U,H,R 

3A0  PRINT  TAB (8) »  Ml #  T1 

380  NEXT  Y 

1000  END 
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"WAKEFL07"  PROGRAM 


1.  Purpose:  To  compute  completely  asymmetric  (P  =  0)  wake  profiles  of 

velocity,  density,  etc.  normalized  so  that  the  maximum  value 
of  each  profile  is  1.0.  This  is  done  for  a  fixed  x',  M& 
and  Tw/Tq.  The  program  requires  inputs  from  the  WAKEFL06  * 

program  in  the  form  of  tlie  extreme  (center)  values  of  velocity, 
temperature  and  density. 

2.  Application  Example 


Compute  the  normalized  profiles  for  the  P  !  0  wake  when  x'  °  0.3,  Mg  *  6, 

Tw/Tq  =  0.2  for  -3  <  y'  <  3  in  steps  of  0.05.  For  this  case,  the 
WAKEFL06  program  gives  the  following  center  values: 

(1  .  u_)  =  0.3705 

ue 

(—-!)  -  1-426 

Te  C  » 

(£_  .  i)  =  -0.5879 
pe 

Enter: 

20  LET  X  =  0.3  • 

22  LET  M  =  6 

24  LET  T  =  0.2 

82  LET  U  =  0.3705 

84  LETT  =1.426  1 

86  LET  R  =  -0.5879 

T 60  FOR  Y  =  -3  TO  3  STEP  0.06 
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"WA^ap?"  LISTING 

10  REM  WAKEFL07 

15  DIM  Y<100>  tZ<100>»B<100>»Yl<100>'U<2L00> 

16  DIM  Z 1  <1 0  0  )  ,  D  <  1  0  0  )  ,  Y  2  <  .1 0  0  >  ,  Y  3  <  1 1)  f)  > ,  H  <  1 0  0  > ,  R  <  1 0  0  > ,  M  <  1 0  0  >  ,  T 1  <  1 0  0  > 

20  LET  Xs ♦ 3 

22  LET  M“6 
24  LET  T~ .  2 

26  LET  B1 ♦  2#  <  MA2  >  #  <  T  *-l )  +T-1 
28  LET  C!~~  ♦  2*  ( MA2 ) 

30  LET  P«. 327591 
40  LET  Al“ -  251829592* 

50  LET  A2:=-«  284496736000  0001* 

60  LET  A3-1  ,421413741# 

70  LET  A4~-l ♦ 453152027# 

80  LET  A5~l ♦ 061405429# 

82  LET  U9= .3705 
84  LET  T9~ 1,426 
86  LET  R9=;-  «  5879 

90  PRINT  "NORMALIZED  PROFILES  FOR  COMPLETELY  ASYMMETRIC  <P=0>  CASE" 
110  PRINT 

.12  0  PRINT  "FOR  X'*"  JX 

122  PRINT  "AND  MACH  NQ.="JM 

124  PRINT  "TW/T0==";T 

125  PRINT  "CENTER  1--IJ/UE  IS="iU9 

126  PRINT  "CENTER  T/TE~1  =  "  ,*T9 

127  PRINT  "CENTER  R/RE~1=="  ,‘R9 

128  PRINT 
.129  PRINT 

130  PRINT  "Y',V'Y'/2SGR(X/  >"♦"<  J -U/UE  >*","<  T/TE-1  >  #'V' (R/RE-1  >  *" 

150  PRINT 

.160  FOR  Y 3  TO  3  STEP  ♦  06 
.170  LET  Z=SGR<X)+Y/<2*SGR<X>  ) 

175  IF  Z<0  GOTO  192 
180  LET  B~1 / ( 1 +P*Z ) 

190  LET  Y1 --  ( A1XB+A2*  <  BA2 )  +A3*  <E!A3>+A4#<E:A4>+A5*(BA5))  #EXP (- <  ZA2  >  > 

191  GOTO  200 

.192  LET  Z=-Y/(2*SQR<X)  )-SOR(X) 

194  LET  B“l/( 1+P*Z ) 

196  LET  Yl~2~ ( A1#B+A2*  <  BA2 )+A3* ( BA3 ) +A4# ( BA4  >  +A5# ( BA5 ) ) *EXP <  - (ZA2) ) 

20  0  LET  U  ~  «  5#EXP ( X+Y ) #Y1 

2.10  LET  Z1-2#SQR < X ) +Y/ < 2*SQR < X )  ) 

220  IF  Z1<0  GOTO  260 
230  LET  D*1/<1+P*Z1> 

240  LET  Y2~(A1#D+A2#<DA2>+A3*<DA3)+A4*<DA4>+A5*<DA5) )#EXP(- <Z1A2> ) 
250  GOTO  290 

260  LET  Z1=-Y/<2*SGR<X> )-2*SOR<X> 

270  LET  D*1/<1+P*Z1> 

280  LET  Y2~2- ( A.l*D+A2* <DA2)+A3#<DA3)+A4#(DA4)  +A5* < DA5  )  )*EXP ( -  (Zl  A2 >  > 
290  LET  H~,5*C*EXP<4#X+2*Y)*Y2+B1*U 

300  LET  Y3-Y/(2*SGR<X) ) 

310  LET  Rb-H/(H+1> 

320  LET  Ml~ < 1~U ) #M*SGR < 1+R ) 

330  LET  Tl" <  < H+l ) * ( 1  +  ♦ 2* ( Ml A2 > ) / < 1  +  ♦ 2* < MA2 ) > ) - 1 

350  PRINT  Y,Y3,U/U9,H/T9,R/R9 
380  NEXT  Y 
390  PRINT 

400  PRINT  "#  NORMALIZED  WITH  THE  APPROPRIATE  CENTER  VALUE" 

1000  END  _ 
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"WAKEFL07"  OUTPUT 


C0MPLETfc~1-y  ASYMMETRIC  <F'=0)  case 


FOR  X'*  .3 
AND  MACH  NO , ~  6 
TW/TO*  ,2 

CENTER  1-iJ/UE  IS*  .3705 
CENTER  T/TE-1*  1.426 
CENTER  R/RE-1--.5879 


Y* 


-3 

-2 

.94 

-2 

.88 

-2, 

.82 

““  £  - 

.76 

-2. 

.7 

-2, 

64 

-2. 

58 

-2. 

52 

r> 

X,  « 

46 

Y'/ZSQRt  X '  > 

a~u/UE>* 

(T/TE-1)* 

(R/RE-l)* 

-2 . 73861 

-2 , 68384 

-2.62907 

-2.S743 

-2.51952 

-2.46475 

-2.40998 

-2.35521 

-2.30044 

-2.24566 

♦381215 

.192365 

♦204186 

.216714 

.229983 

.244032 

.258895 

.27461 

. 291209 

.308726 

. 327995 

.34557 

.363893 

. 382967 

.402795 

.423371 

.444685 

.466722 

.489456 

.512857 

.542051 
.561508 
.581108 
.60081 
.620567 
.640334 
.  660 061 
.6797 
.6992 
.7185033 
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